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ABSTRACT 

We report high spatial and spectral resolution measurements of masers towards R Aqr 
and HI— 36, both of which are examples of the sub-class of symbiotic stars that contain 
a long-period Mira-type variable. Our observations have resulted in the first detection 
of 321-GHz H2O maser action towards a symbiotic Mira — R Aqr. Comparison with 
simultaneous 22-GHz H2O maser data suggests the masers do not have the same 
properties as those in the circumstellar envelopes of field Miras. R Aqr's 22-/321-GHz 
peak flux density and luminosity ratios are low, as is the line-width ratio. Continuum 
and spectral-line maps indicate that the 22-GHz maser and free-free emission are 
aligned. Three mechanisms can reproduce the data with varying degrees of success. 
All three lead naturally to normal levels of maser emission in SiO and 321-GHz H2O 
and anomalously weak OH and 22-GHz H2O masers. In the most convincing model, 
UV radiation and a fast wind from the companion remove the Mira's envelope of dusty, 
molecular gas, leaving a relatively small cavity of dense, neutral material within a large, 
ionized nebula. Excitation temperatures suggest that 321-GHz masers are normally 
excited close to the Mira whilst 22-GHz masers are more remote; in R Aqr, therefore, 
the 22-GHz masers do not form under optimum conditions. Instead, we see weak and 
narrow lines that form closer to the Mira, consistent with our high-resolution maps. 
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1 INTRODUCTION 

R Aqr is the closest known symbiotic Mira. It is a binary 
system in which a late-type, long-period variable giant loses 
matter to a post-AGB dwarf. Existing models have failed 
to account convincingly for many of the simplest features of 
the R Aqr system: the binary period, the method of mass 
transfer, the nature of the hot, companion star and its out- 
burst mechanism. Still less is known about the more unusual 
phenomenon on display: the transient jet-like prominence 
(Wallerstein & Greenstein 1980; Sopka et al. 1982) and the 
arcmin-scale nebulae (Solf & Ulrich 1985). 

Various studies have shown that the molecular emis- 
sion lines of symbiotic stars are very different from those 
expected from isolated Miras. For example. Bowers & Ha- 
gen (1984) stated that 75 per cent of Miras in the solar 
neighbourhood have 22-GHz H2O maser emission, whereas 
Seaquist, Ivison & Hall (1995 — SIH95) reported that only 
two of the twenty-five known Mira-type symbiotic stars sup- 
port 22-GHz H2O masers. R Aqr and HI— 36 are the unusual 
duo. The reasons why these systems are peculiar are not well 
understood, but it seems likely that it is because of the rel- 



atively large size of the binary orbit in the case of HI— 36, 
and because of the low luminosity of the hot companion star 
in the case of R Aqr (compared with other symbiotic objects 
in the samples observed by SIH95 and Schwarz et al. 1995). 

R Aqr and HI— 36 also harbour amongst the brightest 
known SiO masers. R Aqr has been detected in the 43-GHz 
J = 1 — lines of both the t; = 1 and v = 2 states (Martinez 
et al. 1988), with SiO fluxes at 43 GHz averaging around 
100 Jy over several pulsational periods, with the fiux and ve- 
locity of individual features varying considerably from one 
period to the next. It has also been detected in the v = 1, 
J = 2-1 transition at 86 GHz (Hollis et al. 1990). Schwarz 
et al. show that the peak of the 86-GHz emission averages 
about 200 Jy, but varies with stellar phase and between dif- 
ferent periods in both flux and velocity (like the 43-GHz 
emission). Finally, R Aqr is also one of a handful of systems 
to display the J = 7 — 6 SiO transition (Gray et al. 1995). 

In a recent survey of twenty-two circumstellar envelope 
22-GHz H2O maser sources, Yates, Cohen and Hills (1995) 
showed that 70 per cent had observable sub-mm H2O emis- 
sion. The velocity extent of 321-GHz maser emission was 
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Table 1. Parameters for four epochs of 22-GHz observations of R Aqr. 



UT Date 


Array 


Synthesized beam 


Dominant peak 


Position (B1950) 




configuration 


size and PA 


velocity (kms^-'-) 


a(23'^41') 


5(-15°33') 


1993 May 16 


CnB 


1.8" X 0.6", 124° 


-26.6 ± 0.2 


14.27 ±0=.02 


43.30 ± 0".15 


1994 May 20 


BnA 


0.2" X 0.1", 55° 


-26.8 ± 0.2 


14.27 ±0^01 


43.35 ± 0".03 


1995 June 8 


D 


4.3" X 2.6", 154° 


-24.7 it 0.2 


14.26 ±0=. 03 


43.19 ±0". 60 



less than that at 22-GHz and the 22-GHz peak flux density 
was typically twice that at 321-GHz. 

Here, we report the first detection of the J = IO29 — Sse 
ortho 321-GHz H2O maser transition towards a symbiotic 
Mira system. Poor observing conditions prevented any at- 
tempt to detect the para J = 615 — 422 325-GHz H2O 
maser transition which sits in a deep atmospheric absorption 
trough. Our sub-mm observations were performed using the 
James Clerk Majcwell Telescope, whilst observations of the 
J = 616 — 523 22-GHz H2O maser transition were obtained 
contemporaneously using the Very Large Array. 



2 OBSERVATIONS 

2.1 James Clerk Maxwell Telescope 

The JGMT observations were carried out at two epochs. 
A liquid-helium-cooled, single-channel SIS mixer receiver, 
B3i, was used, with a digital autocorrelation spectrometer 
backend providing a bandwidth of 250 MHz and a channel 
spacing of 0.156 MHz. The rest frequency adopted for the 
H2O J = IO29 - 936 transition was 321.22564 GHz, with 
velocities of —25.0 and — 125. km s"'^ relative to the LSR 
for R Aqr and HI— 36, respectively. 

During 1994 May 01, observations were hampered by 
weather conditions unsuited to observing near the 325-GHz 
atmospheric water absorption feature — the opacity at 
230 GHz was ~ 0.1, giving a 321-GHz system temperature 
of around 1500 k. We integrated for Ihr, achieving a 3-cr 
upper limit on the integrated line flux of 13.8 Jy kms"'^. 

At the second epoch, 1995 June 05, Tgys was ~ 1000 k. 
After integrating for 5.25 hr (exclusive of overheads such as 
the time spent nodding the telescope to alternate the signal 
and reference beams), we obtained a detection of the 321- 
GHz transition towaxds R Aqr. 
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Figure 1. The 22-GHz spectra observed towards R Aqr using 

the VLA during 1993 May (lower, from Ivison ct al. 1994), 1994 
May (middle) and 1995 June (upper). Pulsational phases, where 
4> = 0.00 corresponds to maximum visual light, arc (p = 0.33, 0.28 
and 0.27 (Chinarova, Andronov & Schweitzer 1996). Improved 
velocity resolution is apparent in the upper panels. 



2.2 Very Large Array 

Including data reported by Ivison, Seaquist & Hall (1994), 
interferometric observations at 22 GHz have been carried out 
using the VLA on four occasions. Here, we describe the latter 
three epochs. 

The telescope was in the A configuration during 1994 
May 01-02 when data were obtained for HI— 36 and R Aqr 
covering the 1612- and 1667-MHz OH lines simultaneously 
with two IF pairs (for HI— 36: 195-kHz bandwidth, with 
0.28-kms~^ channel separation, for 150 min; for R Aqr: 781- 
kHz bandwidth, with 1.13-kms^^ channel separation, for 
120 min). Data was also obtained for HI— 36 covering the 
22-GHz H2O line (6.25-MHz bandwidth, with 1.32-kms"^ 
channel separation, for 45 min). The HI— 36 data were phase 



referenced against 1748—253 and 1741—312 at the low and 
high frequencies, respectively, whilst the R Aqr data were 
referenced against 2345—167. 

During 1994 May 20, the VLA was used in BnA config- 
uration (i.e. B configuration, with some long N-S baselines) . 
Data were obtained for R Aqr covering the 22-GHz H2O 
line (3.125-MHz bandwidth, with 0.33-kms~^ channel sep- 
aration, for 140 min). A 22-GHz continuum map of R Aqr, 
with 100-MHz bandwidth, was also obtained. 

Our final set of VLA observations were performed, using 
the D configuration, during 1995 June 08 and were contem- 
poraneous with the JGMT measurements described earlier. 
Data were obtained covering the 22-GHz H2O line for both 
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Figure 2. BnA-configuration continuum map of R Aqr at 
22 GHz. The map is marked with the position of the 22-GHz H2O 
maser and a representation of the restoring beam. The peak flux 
density is 12.5 mjybcam"^. Contours arc plotted at integer val- 
ues of the 3-cr level (0.85 mjy bcam^-'^). Sec Ivison et al. (1995a) 
for a line-free map extracted from the 1665- and 1667-MHz OH 
data cubes with the position of the 22-GHz H2O maser marked. 
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Figure 3. The 22-GHz (upper) and 321-GHz (lower) spectra ob- 
served towards R Aqr during 1995 June (0 = 0.27). The 321-GHz 
data from JCMT have been fitted with a Gaussian profile and 
are shown after subtraction of a 2nd-order polynomial fit to the 
baseline. 



HI— 36 and R Aqr (1 hr each, after overheads) with the same 
bandwidths and channel separations as used in 1994 May. 

For all the VLA observations, absolute flux calibration 
was performed using 0134-1-329 and 1328-1-307. 



3 OBSERVATIONS OF R AQUARII 
3.1 The H2O masers 

Fig. 1 shows all the observations of R Aqr obtained so far 
with the VLA at 22 GHz. During our initial maser-line sur- 
vey (reported by Ivison et al. 1994; SIH95), a relatively poor 
velocity resolution was employed, both to improve our sen- 
sitivity to weak lines and to maximise velocity coverage. At 
subsequent epochs we used narrower channels, which is ap- 
parent from the improved velocity resolution of spectra in 
the upper panels of Fig. 1. 

The 22-GHz H2O line from R Aqr shows temporal vari- 
ations similar to those seen in its SiO maser lines (Gray 
et al. 1998). The full velocity width of the 22-GHz H2O 
line (8kms~^) is also the same as that of the 86-GHz 
D = 1,J = 2 — 1 SiO line observed by Schwarz et al. (1995), 
confirming the speculation of Ivison ct al. (1994). 

In terms of its position, the H2O maser is coincident 
with the centroid of the 22-GHz continuum emission (Fig. 2) 
(< 50 milliarcscconds, or < 11 au for a distance of 220 pc — 
van Belle ct al. 1996). These data were obtained contem- 
poraneously. During 1994 May, the full velocity extent of 
8 km s~^ and the presence of several spectral features sug- 
gests that the 22-GHz maser region of R Aqr is complex 
and made up of several discrete masing clumps. The 0.22 x 
0.11 mas^ restoring beam did not resolve any structure, but 
it does allow us to confirm that the emission mechanism is 
definitely non-thermal, since Tb > 12,000 K. 



At the other two epochs, there was one dominant peak. 
The velocities of the dominant line peaks during 1993 and 
1994 were separated by only 0.2 kms^^, but a 2.1-kms~^ 
redward shift is evident in the 1995 data. The full profile 
width is similar at all three epochs, always covering the —22 
to — 30kms~^ v\sT interval. 

Fig. 3 shows the 22- and 321-GHz spectra obtained dur- 
ing 1995 June for R Aqr. The 321-GHz spectrum is noisy, so 
it would be a mistake to draw too many conclusions from the 
line profile or width. The best Gaussian fit to the data yields 
a FWHM line width of 8 ± 2kms~'^ with a central velocity 
of —26 ± Ikms"^ (slightly blueward and broader than the 
dominant 22-GHz line peak — see below). The integrated 
fiux density is 10.5 ± 2.5 Jykms"^ and the corresponding 
photon luminosity is (3.06 ± 0.73) x lO''^ s"\ 

The upper plot of Fig. 3 shows the contemporanous 
22-GHz spectrum. As with all our VLA data, the line emis- 
sion stretches from —22 to — 30kms~^; the dominant peak 
(0.37 Jy after subtracting off the free-free continuum con- 
tribution) is at visT ~ — 24.7kms~^; its FWHM is around 
1.8 ± 0.1 kms^^ whilst the integrated flux density of the en- 
tire profile, including the pedestal, is 0.73 ±0.05Jykms~^ 
or (6.4 ±0.4) X 10*° photons"\ 

3.2 Comparision with the H2O masers of field 
Miras 

The recent survey of the profile characteristics of 22- and 
321-GHz maser emission from twenty field Miras (Yates et 
al. 1995) gave a mean ratio of 2.2 for the peak flux densities 
of the 22- and 321-GHz lines and a mean luminosity ratio of 
1.5. There was no apparent correlation between these ratios 
and the mass-loss rate or the IRAS colours, though the line 
luminosities did correlate with both of these parameters. 
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For R Aqr we measured a peak flux density ratio of 
0.24 ± 0.05 and a photon luminosity ratio of 0.21 ± 0.05. 
These he at the lowest bound of the values observed by 
Yates et al. (1995). The only objects in the Yates et al. 
sample that rival those ratios arc o Ccti, U Her, R Leo and 
R Cas; the former is mildly symbiotic - a wind-accreting 
companion star is separated from the Mira by around 80 au 
(Karovska, Nisenson & Belefic 1993); all four have relatively 
low mass-loss rates. R Aqr also has a very low mass-loss rate 
(0.06 X lO^^Moyr^S Soaquist, Kroguloc & Taylor 1993). 

As was mentioned earlier, Yates et al. (1995) found 
no correlation between the 22-/32 1-GHz maser-line pho- 
ton luminosity ratio and the mass-loss rate for a sample 
of 13. We have re-analysed the data presented by Yates 
et al., supplementing them with the result for R Aqr and 
treating upper and lower limits for the 22-/321-GHz maser- 
line photon luminosity ratio as though they were detec- 
tions (the effect of which is to underestimate the degree of 
correlation, since the systems with low mass-loss rates are 
those with upper limits and the systems with high mass- 
loss rates are those with lower limits). With a sample of 
19, the resulting correlation coefficient is 0.66, indicating a 
99.8 per cent probability that the 22-/321-GHz maser-line 
ratio varies systematically with the mass-loss rate such that 
logio(L22/L32i) = 0.32 0.88 logio(M/10-^ M© yr'^). This 
would be expected if the maser transitions are inverted at 
different distances from the central star and a similar trend 
was tentatively reported by Menten & Melnick (1991). 



3.3 The OH lines: upper limits 

1612- and 1667-MHz spectra of R Aqr were obtained using 
the VLA during 1994 June. The 2-hr integration resulted in 
a substantially lower noise level than that achieved by Ivison 
et al. (1994); still, neither of the OH lines were detected. The 
noise level (in a region covering TJisr ~ —96 to +46kms~^) 
was 2.7mJy for both spectra (with l.l-kms~^ channels). 

If the line tentatively identified by Ivison et al. (1994) at 
— 14kms~^ were a real feature (and OH masers are usually 
variable on much larger timescales than masers that form 
closer to the LPV) then we would have detected it at the 
~ 30-cr level. Our limit on the integrated flux (assuming 
a line width of 8.6kms"^) is S{3a) < 0.025Jykms"\ or 
< 1.38 X 10"^" Wm-^ and we conclude that the feature 
described by Ivison et al. (1994) was spurious. 



4 OBSERVATIONS OF Hl-36 

Fig. 4 shows the OH and H2O spectra of HI— 36 obtained 
during 1994 May. Both the spatial and the velocity resolu- 
tion were better than those obtained by Ivison et al. (1994) 
and SIH95 during the survey-style observations, allowing us 
to rule out a thermal origin for the OH and H2O line emis- 
sion (Tb > 9 X 10^ and 3 x 10"* K, respectively). Despite this, 
the longer integration times reported here meant that the 
noise values in the OH and H2O spectra were several times 
lower than in our 1993 spectra (see Table 2). 

The 1612-MHz OH line centroid {vi^r = -124.5 ± 
l.Okms"^) is displaced by — 5.5kms~^ from the mainline 
centroid, though a faint red wing stretches to the centre of 
the mainline. The full widths are ~ 8kms~^ for both of the 
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Figure 4. Masers observed towards HI— 36 using the VLA in A 
configuration during 1994 May 01. 
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Figure 5. A-configuration 8.4-GHz isophotal map of HI— 36 ob- 
tained during 1991 Sept). The positions of the OH and 22-GHz 
H2O masers are marked and the size of the restoring beam is 
shown. Contours are plotted at 1, 2, 4 ... 128 times the 3-a level 
of 0.25 mjy beam~^ . 
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Table 2. Parameters for the VLA A-configuration observations of HI— 36 during 1994 May 01. 



5 





Line 


Velocity 


Spectral 


Beamsize and 


Central 


Integrated flux 






resolution 


noise (mjy) 


PA (" and °) 


velocity 


density (Jykms^-'-) 


OH - 


1612 MHz 


0.28kms"l 


5.2 


2.2 X 0.6, 172 


-124.5 ± l.Okms-l 


0.71 ±0.05 


OH - 


1667 MHz 


0.27kms-l 


4.5 


2.1 X 0.6,172 


-119 ± l.Okms-i 


0.19 ±0.03 


H2O 


- 22 GHz 


1.32kms-l 


4.0 


0.29 X 0.06, 22 


-121 ± 1.5kms-i 


0.16 ±0.04 



lines, whilst the FWHM are ~ 4 and ~ 6 kms"^ at 1612 and 
1667 MHz, respectively. The velocity of the strongest peak 
at 1612 MHz corresponds with the 1993 May and July val- 
ues (to within the uncertainties, Ivison et al. 1994). This is 
also true of the H2O Une (which is at only 25 per cent of its 
1993 intensity). 

If the observed 1612-MHz peak is the blueward compo- 
nent of the usual OH pair, then we can estimate a termi- 
nal wind velocity of lOkms"^ from its displacement from 
Hl-36's SiO and H92a lines (Allen et al. 1989; Bastian 
1992); furthermore, if the unseen redward peak has been ob- 
scured by optically thick, ionized gas (an idea first proposed 
by Allen et al. 1989 — see §5.2.3) then t > 3.6 at 1612 MHz, 
based on the spectral rms redward of the observed line. 



5 ON THE R AQR SYSTEM 

5.1 Radial velocity measurements in the radio 

Several estimates of R Aqr's orbital period have been pro- 
posed; none, however, have been robust. Table 3 lists the 
orbital periodicities that have been suggested to date, whilst 
more recent radial velocity (RV) data have been presented 
by Wallerstein (1986) and Hinkle ct al. (1989). The RVs do 
not strongly support any of the suggested periodicities. Hin- 
kle et al. state that the period must be between 18 and 21 yr 
or longer than 32 yr. Most solutions for periods of over 32 yr 
require e > 0.3. 

Monitoring the position and velocity of the R Aqr wa- 
ter maser could reveal the system's elusive binary period. 
The VLA offers both high velocity resolution and high po- 
sitional accuracy with which to track a Mira. Table 1 shows 
the measurements from our first three epochs of monitor- 
ing; so far, there are no appreciable positional shifts and 
only one appreciable velocity shift. The problem could be 
that the timescale for the rise and fall of individual maser 
features (which is weatly correlated with the pulsational pe- 
riod) is shorter than the orbital period. This is compounded 
by strong, new maser features (such as that seen at the 0.4- 
Jy level in 1995 June) which can obscure weaker features 
such as the 0.1- Jy peak of 1994 May (see Fig. 1) . 

The 8-kms~^ pedestal is the most enduring feature of 
the maser spectrum, and its stability at the Avisr < 1 kms~^ 
level over a period of two years implies a low orbital incli- 
nation, or a high Mgiant : Afdwarf ratio, or an orbital period 
that is measured in decades, which is consistent with opti- 
cal data. If the pedestal is the result of thermal emission (at 
the 10-mJy level) from near the photosphere of the rotating 
LPV then its FWZM translates into a lower limit for the ro- 
tational period. The period will approach the lower limit if 
intrinsic line broadening is negligable, i.e. if the width is due 
entirely to Doppler broadening. The maximum photospheric 



Table 3. Orbital periods suggested in the literature. 
Periodicity Description 

/yr 

26.7 RV measurements (Merrill 1950). Later refuted 

by Jacobson & Wallerstein (1975). 
24.7 Same RV data as Merrill (1950), used by 

Kurochkin (1976). 
44 Prom two periods of peculiaj: activity in the 

Mira's LC that are suggestive of its eclipse 

(Willson et al. 1981). 



line-of-sight velocity implied by the pedestal is 4kms~ , 
well within the scatter of the observed rotational velocities 
for late- type giants (vsini = 1 — 9kms~^ for K giants, De 
Medeiros, Melo & Mayor 1996). This yields a rotational pe- 
riod of > 17yr for an LPV radius of 500 ± 100 R© (Haniff, 
Scholz & Tuthill 1995; van Belle et al. 1996). 

5.2 Mechanisms to explain the simultaneous data 

5.2.1 Disruption by the white dwarf secondary? 

Any periodicity in excess of 8 yr would yield a binary sepa- 
ration which would put the hot star in the H2O maser zone 
(10 — 20 au), i.e., the companion star would orbit in the re- 
gion normally occupied by the masing molecular gas. The 
companion could then reduce amplification by inducing tur- 
bulent motion and disrupting the velocity-coherent paths in 
the outer envelope where the bright 22-GHz masers form 
(Hall et al. 1990a, 1990b). The interaction of the accretion 
flow and the Mira's wind could cause large velocity gradients 
and turbulence in the circumstellar envelope, which reduce 
the velocity coherence path (gain length). 

5.2.2 Consistency with SIH95 model? 

S1H95 put forward another model that appears to be consis- 
tent with the results of not only their survey of maser-line 
transitions towards symbiotic Miras, but also with the de- 
tails of individual systems such as R Aqr and HI— 36. In their 
model, UV radiation and a fast wind from the white dwarf 
companion dissociate and sweep away the dusty, molecular 
gas shed by the Mira. Even the molecular material shielded 
by the Mira from direct irradiation becomes dissociated as 
the result of Rayleigh scattering in the Mira's wind. How- 
ever, since the R Aqr hot component has a relatively low 
luminosity (and the nebula is ionisation bounded), the ef- 
fectiveness of direct and indirect molecular dissociation by 
UV irradiation would be much reduced; even on the irrar 
diated hemisphere of the Mira, an extensive neutral cavity 
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may remain beyond the ionization front of the hot dwarf. 
Thus we have an environment where 22-GHz H2O masers 
may survive, in keeping with our observations. 

The energy levels of the 22-GHz and 321-GHz transi- 
tions arc 643 and 1861 K above the ground state, respec- 
tively. Yates, Field & Gray (1997) suggest that 321-GHz 
masers are excited closer to the Mira than 22-GHz masers. 
Given the absence of a suitable dusty, molecular environ- 
ment in the outer reaches of the R Aqr Mira's wind, the 
observed 22-GHz H2O rnascr lines will be weak and narrow 
in comparison to those observed towards field Miras. 

Note also that our maps (Fig. 2) seem to indicate that 
the 22-GHz mascr and the ionized gas arc co-spatial (§3.1). 
This makes little sense unless the molecular gas is protected 
in a neutral cavity, much closer to the LPV than would 
normally be anticipated. The superimposition can be ex- 
plained if the maser is within a cocoon of ionized gas and 
the strongest contribution to the free-free radiation compo- 
nent comes from the densest portion of the ionized wind, 
which is near the ionisation front and close to the Mira. 

5.2.3 Line obscuration at low frequencies? 

This was proposed Allen et al. (1989), and relies on the 
fact that the LPV's wind is optically thick at low frequen- 
cies. The spectral energy distributions of HI— 36 and R Aqr 
show the turnover to optically thin free-free emission is at 
~ 10 GHz (Ivison, Hughes & Bode 1992; Ivison et al. 1995b). 
Below 10 GHz, we can expect the line emission to be hidden 
from us by the optically thick, ionized gas. 

The rnascr data support Allen's idea: first, the absence 
of OH maser emission towards R Aqr (see §3.3); second, only 
one component of the normal 1612-MHz OH mainline pair 
is seen towards HI— 36, which suggests a large optical depth 
between the front and back of the masing shell in the line of 
sight through the system. Third, the observed lines become 
progressively more intense with increasing line frequency, 
suggesting lower obscuration: the OH line is absent, the 22- 
GHz emission is weak and the line emission at 43, 86, 299, 
301 and 321 GHz is strong. There is no evidence, however, 
that the opacity at 22 GHz is sufficient to move the 22-/321- 
GHz line flux and luminosity ratios closer to those typically 
observed towards field Miras (a 22- /321-GHz line luminosity 
ratio of 1.5 would require r ~ 2 at 22 GHz). 



6 CONCLUDING REMARKS 

We have presented maser-line spectra and maps of the sym- 
biotic Miras, R Aqr and HI— 36, together with contempora- 
neous, high-resolution continuum maps. Our measurements 
have resulted in the first detection of sub-mm water maser 
emission towards a symbiotic system (R Aqr) and we have 
discussed three mechanisms that could cause the observed 
line properties in the geometries implied by the continuum 
and spectral-line maps. 

With the data available to us, it is not possible to dis- 
criminate between three possible scenarios, although one of 
these requires conditions that arc not supported by other 
measurements. The mechanisms that we favour inhibit the 
formation of OH and 22-GHz H2O masers in their natural 
environments, employing either the orbital motion of the 



companion or its radiation field and fast wind. The third, 
which reduces maser-line intensities via the free- free opacity 
of circumstellar, ionized gas (and thus has an effect that is 
dependent on the rest frequency of the meiser transition) is 
somewhat less credulous since to reduce the intensity of the 
22-GHz maser line by at least an order of magnitude rela- 
tive to the 321-GHz line would require the free-free emission 
to remain optically thick at that frequency, which does not 
appear to be the case (see Ivison et al. 1992, 1995b). 
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